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' We investigate the phase diagram of color superconducting quark matter with strange quarks 

Q \ (2SC+S quark matter) in beta equhibrium at zero as well as finite temperatures within a Nambu- 

Jona-Lasinio model. The variational method as used here allows us to investigate simultaneous 
1^-^ ' formation of condensates in quark-antiquark as well as in diquark channels. Color and electric 

charge neutrality conditions are imposed in the calculation of the thermodynamic potential. Medium 
dependance of strange quark mass plays a sensitve role in maintaining charge neutrality conditions. 
At zero temperature the system goes from gapless phase to usual BCS phase through an intermediate 
^ ' normal phase as density is increased. The gapless modes show a smooth behaviour with respect 

CO , to temperature vanishing above a critical temperature which is larger than the BCS transition 

temperature. We observe a sharp transition from gapless superconducting phase to the BCS phase 
04 ' as density is increased for the color neutral matter at zero temperature. As temperature is increased 

, this however becomes a smooth transition. 

T— I . 

■>!::;j- . PACS numbers: 12.38.Mh, 24.85.+p 

o 

p • I. INTRODUCTION 



X 



Color superconductivity has become a compelling topic in QCD during the last few years. At sufficiently high 
baryon densities, when nucleons get converted to quark matter, the resulting quark matter is in one kind or the other 
of the many different possible color superconducting phases at low enough temperatures T| . The rich phase structure 
is essentially due to the fact that the quark quark interaction is not only strong and attractive in many channels but 
■ also many degrees of freedom are possible for quarks like color, flavor and spin so that various kinds of BCS pairing 
are possible. Studying properties of color superconducting phases in heavy ion collision experiments seems unlikely 
in the present accelerators as one cannot avoid producing large entropy per baryon in heavy ion collisions and hence 
cannot produce the dense and cold environment that is needed to supoort the formation of superconducting phases. 
However, in the future acclerator facility planned at GSI for compressed baryonic matter experiments, one possibly 
can hope for observing fluctuations signifying precursory phenomena of color superconducting phase [2l ■ 

On the other hand, it is natural to expect some color suprconducting phase to exist in the core of compact stars 
where the densities are above nuclear matter densities and temperatures are of the order of tens of keV. However, 
to consider quark matter for neutron stars, color and electrical charge neutrality conditions need to be imposed for 
the bulk quark matter. Such an attempt has been made in Ref.i as well as in Ref. where the lighter up and 
down quarks form the two flavor color superconducting (2SC) matter while the strange quarks do not participate 
in pairing. A model independent analysis was done in Ref. |3] that is valid in the limit rus « fi and A ^ m^^/ ji, 
where, A is the pairing gap and /i is the quark chemical potential. It has been shown, based upon comparison of free 
energies that such a two flavor color superconducting phase would be absent in the core of neutron stars 3] . Within 
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Nambu Jona-Lasinio (NJL) model in Ref. [J| it has been argued that such conclusions arc consistent except for a 
small window in density ranee where superconducting phase is possible. There have also been studies to include the 
possibility of mixed phases Q of superconducting matter demanding neutral matter on the average. Recently it was 
observed that imposition of neutrality conditions lead to pairing of quarks with different fermi momenta giving rise 
to gapless modes Within a Nambu JonaLasinio model the two flavor superconducting quark matter (2SC) was 

shown to exhibit gapless modes (g2SC) arising due to difference in the fermi momenta of the pairing quarks when 
charge and color neutrality conditions are imposed. Superconducting quark matter with strange quarks (2SC+s) was 
shown to exhibit gapless superconductivity (g2SC+s) within a window of about 80 MeV in baryon chemical potential 
0. The number densities of the strange quarks were seen to be as large as about 40% of the number density of e.g. 
the up quarks in the gapless phase 7]. Temperature effects on the gapless modes were also studied for two flavor 
quark matter in Ref.|3]. 

We had applied a different approach to study the problem in Ref. 0, 0| . We considered a variational approach 
with an explicit assumption for the ground state having both quark-antiquark as well as diquark condensates. The 
actual calculations are carried out for the NJL model with a vector interaction such that the minimisation of the free 
energy density determines which condensate will exist at what density. Charge neutrality conditions were introduced 
through the introduction of appropriate chemical potentials. We note here that possibility of diquark condensates 
along with quark-antiquark condensates has also been considered in Ref. 0,Ei.l2, 13] . 

In the present work we generalise the variational approach of Ref. to study the color superconductivity as well 
as the chiral symmetry breaking to include the effects of temperatures. This will be particularly relevant for the 
physics of proto neutron stars. The gapless modes in the two flavor color superconductivity were investigated at 
finite temperatures Q. However, in the present work we take the analysis further in the sense that the diquark and 
the quark-antiquark condensates are uniquely calculated selfconsistently and elaborate the transition from gapless 
phase and BCS phase as the density is increased for different temperatures. Further, since the realisation of color 
and charge neutrality is nontrivial for physically relevant values of strange quark mass, we retain all orders in rus in 
the calculation of the thermodynamic potential. This of course is important to understand the phase structure of 
quark matter at densities relevant to the interior of neutron stars. The temperature dependence and the nature of 
the gapless modes as a function of the baryon density shall be particularly relevant for cooling of the neutron stars 
through neutrino emission .14] . 

We organize the paper as follows. In the next subsection we discuss the ansatz state with the quark-antiquark as 
well as the diquark condensates. In section ^1 we consider the Nambu Jona-Lasinio model Hamiltonian and calculate 
its expectation value with respect to the ansatz state to compute the thermodynamic potential. We minimise the 
thermodynamic potential to calculate all the ansatz functions and the resulting mass as well as superconducting gap 
equations here. In section Hill we discuss the results of the present investigation. Finally we summarise and conclude 
in section Hvl 



A. An ansatz for the ground state 

To make the notations clear, let us note first the quark field operator expansion in momentum space given as |l5lll6j 



1 



(27r)3/2 



[Uoik)q^{k) + yo(-k)g"(-k)] e^-^-dk, (1) 



where. 



The operators q'^ and (f are the two component particle annihilation and antiparticle creation operators respectively 
which annihilate or create quanta acting upon the perturbative or the chiral vacuum ]0). We have suppressed here 
the color and flavor indices of the quark field operators. The function 0°(k) in the spinors in Eq.(|2l are given as 
cot0° = mi/]k], for free massive fermion fields, i being the flavor index. For massless fields 0°(]k]) — tt/2. 



3 



We now write down the ansatz for the variational state as a squeezed coherent state involving quark antiquark as 
well as diquark condensates as given by ^ 

\n) =Ud\vac) =UdUQ\0), (3) 
Here, Uq and Ud are unitary operators creating quark-antiquark and diquark pairs respectively. Explicitly, 

Uq = exp (^J q"'{k)\(T ■ k)h,{k)q"\k)dk~h.c}j . (4) 

In the above, /ii(k) is a real function of |k| which describes vacuum realignment for chiral symmetry breaking for 
quarks of a given flavor i. We shall take the condensate function /i(k) to be the same for u and d quarks and /i3(k) 
as the chiral condensate function for the s-quark. Clearly, a nontrivial hi (k) shall break chiral symmetry. Summation 
over three colors and three flavors is understood in the exponent of Uq in Eq. Q). Similarly, the unitary operator Ud 
describing diquark condensates is given as 



Ud = exp(Bt - Bd) 
where, i?^ is the pair creation operator as given by 

bI= [ knk)V/'nk)d',(-k)te,,e3,fc + g-(k)r/r(k)e''.(-k)e,,e3aJ dk. 



(5) 



^ (6) 

In the above, i,j are flavor indices, a,b are the color indices and r{= ±1/2) is the spin index. The operators g(k) 
are related to q°(k) through the transformation (/(k) = UQq^(k)UQ^ . As noted earlier wc shall have the quarks of 
colors red and green (a=l,2) and flavors u,d (i=l,2) taking part in diquark condensation . The blue quarks (a=3) 
do not take part in diquark condensation. We have also introduced here (color, flavor dependent) functions P°'(k) 
and /{"^(k) respectively for the diquark and diantiquark channels. As may be noted the state constructed in Eq.jSJ 
is spin singlet and is antisymmetric in both color and flavor. Clearly, by construction /'"(k) = /■'^(k) with i ^ j and 
a ^ h. The ansatz for the ground state has been written down keeping quark-antiquark condensates for the three 
flavors and diquark condensates for two flavors and two colors. 

Finally, to include the effects of temperature and density we next write down the state at finite temperature and 
density |i7(/3,/i)) taking a thermal Bogoliubov transformation over the state |f2) using thermofield dynamics (TFD) 
as described in ref.s [13, Ell- We then have, 



where Ufi^^ is given as 



with 



^Up,^,\^i) =Up^^UdUQ\{)). 



dk. 



(7) 
(8) 
(9) 



In Eq.@ the ansatz functions 0±(k, /3,/i) will be related to quark and antiquark distributions and the underlined 
operators are the operators in the extended Hilbert space associated with thermal doubling in TFD method. In Eq. 
we have suppressed the color and flavor indices on the quarks as well as the functions 9{k, /3, /i). All the functions in 
the ansatz in Eq.jTI) are to be obtained by minimising the thermodynamic potential. This will involve an assumption 
about the effective Hamiltonian. We shall carry out this minimisation in the next section. 



II. EVALUATION OF THERMODYNAMIC POTENTIAL AND GAP EQUATIONS 



Since we shall be dealing with non-asymptotic densities we can not use the weak coupling method. For this reason 
we shall be working in a model in which the interaction among the quarks is simplifled - while still respecting the 
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symmetries of QCD. One natural choice is to model the interaction between the quarks using a four fermion point 
interaction - namely, the Nambu Jona Lasinio model. 
The Hamiltonian is given as 

8 

i,a A=0 

where ■0*'° denotes a quark field of color a and flavor i, and rrii is the current quark mass. As noted earlier, we 
shall assume isospin symmetry with mu=ni^. The interaction term in Ea. pU|l is written down assuming that the 
underlying gluonic degrees of freedom can be frozen into point like effective interactions between the quarks. In 

Ea. (|10|l A"*, A = 1, • • - 8 denote the Gellman matrices acting in the flavor space and A° = y^l/, 1/ as the unit 

matrix in the flavor space. The point interaction produces short distance singularities and to regulate the integrals 
we shall restrict the phase space inside the sphere |p| < A - the ultraviolet cutoff in the NJL model. This form 
of Lagrangian can arise e.g. by Fiertz transformation of a four point current current interaction having quantum 
numbers of single-gluon exchange In that case the diquark coupling Gd is related to the scalar coupling as 

Gd = 0.75Gs- In the study of meson spectrum in this model [20l f2ll | one adds a six point t'Hooft interaction which 
breaks U{1)a symmetry. Although it is straightforward to include such a term, we neglect it here for the sake of 
simplicity. This can have some qualitative consequences as we shall discuss later. 

Using the fact that the variational ansatz state in Eq . arises from successive Bogoliubov transformations one can 
calculate the expectation values of various operators Q . These expressions are used to calculate thermal expectation 
value of the Hamiltonian to compute the thermodynamic potential. A straightforward, but cumbersome calculation 
yields 

e= (H) ^T+Vs + Vd (11) 

arising from the kinetic energy part and the scalar and diquark interaction terms of the Hamiltonian respectively. 
Explicitly, the kinetic energy part in Ea. l|lU|l is given as 

T = (r!(/?,^)|^t(_,c..V + 7V)V|f^(/3,M)) 

2 3 3/. 

= EE / rfk(|k|sin0'+m,cos0,)(l-cos2/i,(k)(l-F''^-i^r)), (12) 

^ i=l a=l 

where, F*" and F^'^ are defined as 

F"(k) ^ sin^ 6'L"(k) + sin^ (Ci'^(k) - sin^ 0'j'{k)) (l - S""^) , (13) 

and, 

Fl'^ik) = sin^ 6i;"(k) + sin^ /j"^ (cos^ 6i;"(k) - 5;"(k)) (l - S""^) . (14) 

Here, we have defined C!? = |e"'e''"'| cos^ 6'*>' and S'^ = le^'e""'] sin^ 6'^'''. The 6"^ term indicates that the third 
color does not take part in diquark condensation. 

The contribution from the scalar interaction term in Ea. HlQ|) turns out to be 



Vs ^ -Gs{n{p,^,)\ ^ [(viA^^)^ - (V^t'a^v)'] \mt^)) - -8G5 E ^h' (i5) 



A=0 1=1,3 

where, 

Il=Yl [ dk cos F'"" - Fi") (16) 



is proportional to (V'V*) ( not summed) and the function (/>,; is defined in terms of the quark-antiquark condensate 
functions as 4>i(k) = 9^° (k) — 2hi{k). 
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In a similar manner, the contribution from the diquark interaction from Ea. (|lU|) to the energy density is given as 

yij = l,2:a,6=l,2 / 

with 

^Wfl '^^^'"''^^^ (^^) (^^^ 

where, S'*"'-'^ has been defined as 

S'*^''^''(k) = sin2r'^(k)cos26l"'^''(k,/3,/i) +sin2/f (k)cos26l!^"'^''(k,/3,^). (19) 

In the above we have defined cos 29^^'-'^ = 1 — sin^ 9!^ — sin^ 9-!^, with i,j = 1,2 being the flavor indices and a,b = 1,2 
being the color indices and i ^ j, a ^ b. 

To calculate the thermodynamic potential we shall have to specify the chemical potentials relevant for the system. 
Here we shall be interested in the form of quark matter that might be present in compact stars older than few 
minutes so that chemical equilibriation under weak interaction is there. The relevant chemical potentials in this case 
then are the baryon chemical potential fis — 3/i, the chemical potential /is associated with electromagnetic charge 
Q = diagi^l"^^ —1/3, —1/3) in flavor space, and the two color electrostatic chemical potentials /is and /ig corresponding 
to C/(l)3 X [/(1)8 subgroup of the color gauge symmetry generated by cartan subalgebra Q3 = diag{l/2, —1/2, 0) and 
Qs = diag{l/Z, 1/3,-2/3) in the color space. Thus the chemical potential is a diagonal matrix in color and flavor 
space, and is given by 

A*ii,a6 — if^Sij + QijfJ.E)Sab + (Qsafc + QsafcAis)^- (20) 

Here, i,j are flavor indices and a, b are color indices. 

The total thermodynamic potential, including the contribution from the electrons, is then given by 

n = T + Vs + VD-{f^N) -^s + n, (21) 

p 

where, we have denoted 

{^^N) = (V^''V.,,a6V'^'') = 2 J2 M^'^C (22) 

i,a 

with fi^°' being the chemical potential for the quark of flavor i and color a, which can be expressed in terms of the 
chemical potentials /i, /i^;, /is and /is using Ea. (|20|l . Further 

^ (2^)^ / '^^^^^'^ ~ ^^^^ 

is proportional to the number density of quarks of given color and flavor. The thermodynamic potential for electrons 
is given as 



l + 27r% 



where we have taken the electron mass to be zero which suffices for the system we are considering. 
Finally, for the entropy density for the quarks we have |17| 



- — ^ V / dk f sin^ 61" In sin^ 6*^" + cos^ In cos^ 9'" 

^ ^ 'i , a 

sin^ 9'^ In sin^ 9!^ + cos^ 9'^ In cos^ 9'^^ . (25) 
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Now functional minimisation the thermodynamic potential f2 with respect to the chiral condensate function /ii(k) 
leads to 

cot0.(k) = l^l (26) 

where, Af^ = nii + SGgll- Substituting this back in Ea. (|16|l yields the mass gap equation as 



where, ej — y k^ + M| being the energy of the constituent quarks of j-th flavor. Clearly, the above includes the effect 

of diquark condensates as well as temperature and density through the functions F and Fi given in Eq.s (|13|) and 
((njl respectively. 

As noted earlier, there are only two independent diquark condensate functions /^^ and /^^ corrsponding to the 
indices u-red and u-green respectively. This is due to the fact that the function with indices d, green is the same as 
that with indices u,red and the condensate function with indices d-red is the same as that with indices u-green. 

Minimisation of the thermodynamic potential fl with respect to the diquark condensate functions /^^(k) and /i2(k) 
yields 

tan2/"(k) = ^— cos(^l^) (28) 
e - Mil 2 

and 



A 



tan2/-(k) = ^^cos(^^^-^) (29) 
e — Hi2 ^ 

where, we have defined A = 4:Gd{I^^'^^ + I^^'^^) with /'"■J'' as defined in Eq.l^Hll- Further, in the above e = (ei+e2)/2, 
Mil = (mii+M22)/2, Mi2 = (/ii2+M2i)/2. It is thus seen that the diquark condensate functions depend upon the average 
energy and the average chemical potential of the quarks that condense. We also note here that the diquark condensate 
functions depends upon the masses of the two quarks which condense through the function cos ((0i — 02)/2). The 
function cos = Mi/ct, can be different for u,d quarks, when charge neutrality condition is imposed. Such a 
normalisation factor is always there when the condensing fermions have different masses as has been noted in Ref. 
[2^ in the context of CFL phase. 
In an identical manner the di-antiquark condensate functions are calculated to be 

tan2A"(k) ^ ,^cos(^l^) (30) 
e + Mii 2 

tan2/-(k) = ^cos(^l^). (31) 

Substituting the solutions for the condensate functions in the expressions for /'"J'' in Ea. pSfl we have the gap equation 
for A given as 



A = 4Gz5(/"-22+/12.21) 

/ oal2,21 r,zill.22 oal2-21 ozill-22 \ 

/ cos 26' cos2t^ cos 26', cos 26', \ ■, , v-i v-/ i , ^ 

rfk A + ^-.n + -.it + ~.it cos^ ^^^-^ (32) 



(27r)3 



In the above, lu'^ = y^A2 cos2((0i - 02)/2) + C±ia = e ± /2*° ,and, cos26'±'^'' has been defined after Eq.lO. 
Finally, the minimisation of the thermodynamic potential with respect to the thermal functions 0±(k) gives 

sin^ e'f = —5-^^ (33) 

^ exp(/3a;!f ) + 1 ^ ' 
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Various w'^'s are given explicitly as follows. 



± s:. 



,21 



-,12 



SeTS 



12 



/ >22 



(34) 



i3 



Also for the strange quark (i=3) u'^ 



g3_|_^3a^ with 



and, finally, for the noncondensing colors oj^ 
£3 — ■\/k2 + M|. As already mentioned, the first index refers to flavor and the second index refers to color. Here 
6e = (ei — £2)72 is half the energy difference of the two quarks which condense and e.g. 5j^^ = (/zn — /Lt22)/2, is 
half the difference of the chemical potentials of the two quarks which condense. Further ui't^ are defined earlier after 
Ea. H32|l . Note that in the absence of imposing the charge neutrality condition all the four quasi particles will have the 
same energy It is clear from the dispersion relations given in Ea. (|34|) that it is possible to have zero modes, i.e., 
— depending upon the values of 6^ and 6^. So, although we shall have nonzero order parameter A, there will 
be fermionic zero modes or the gapless superconducting phase [23!, '2^. We shall discuss more about it section UTTl 

Next, let us focus our attention for the specific case of superconducting phase and the chemical potential associated 
with it. First let us note that the diquark condensate functions depend upon the average of the chemical potentials 
of the quarks that condense. Since this is independent of we can choose to be zero. In that case, = 
/! + 2/3/i_E + /^s/VS = and also /z^^ = fi — l/3fiE + fJ-s/ — fJ-^^ ■ Thus the chemical potentials of the fight fiavors 
become degenerate for both the colors that take part in condensation. This also means that the average chemical 
potential of both the condensing quarks are the same i.e. /2^^ — p}^ — p. = fi + 1 /G^ie + Ms/n/S- For the same reason 
we also have S^^^ = ^_e/2 = d^^^ = 5^. 

With this condition, it is also clear from Eq (|34(l that the quasi particle energies for each flavor becomes degenerate 
for both the colors which take part in condensation. Thus the quasi particle energies now become W-i = ui- + — 
and uj^2 — ^- 



5^, for u and d quarks respectively. Similarly for the antiparticles the energies are given as 



6fj, and uj+2 



Sfj_, for u and d quarks. The gap equation given by Ea. l32|l then reduces to 



A 



8Gd 
(2^ 



dkA 



sin^ ^2) 



UJ+ 



„2 fll 



,2 a2 - 



COS^(- 



(35) 



where, lu± = y/A^ cos'^{{(j)i - 02)/2) + (e± fi)^. 

Now using these dispersion relations, the mass gap equation Ea. H27|) and the superconducting gap equation Ea. (|35|) . 
the thermodynamic potential given in Ea. ()21|l becomes 



^ — ^ud + + ^e, 

where, flud is the contribution from the light quarks and is given as 



(36) 



(27r)3 
4 



A2 



4G 



D 



1 



(CJ_ + UJ+) 



^3 / rfk[log(l + exp(-/3w_,)+log(l + exp(-/3cj+,)] 

2=1,2'^ 



E 



8Gs 



^|dk[v/k^-i(.,+£2) 

1 dk [log(l + exp(-/3(e, - m'^)) + log(l + exp(-/?(e, + /Z^))] 



(37) 
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Similarly the contribution from the strange quarks to the thermodynamic potential, Qs is given as 
6 



(2^) 



^ ^ /'dk[log(l + exp(-/3(e3-M''"))+log(l + exp(-/3(e3 + M''"))] ' 

(38) 



/3(27r)3 ^ y ■ -^v ^-v- ^ ■ ^-'^v- ■ -^v ^-v- ■ ^ 'Ji • gQ 



Finally, the contribution of the electron to the total thermodynamic potentail fie is as given in Ea. l24|l . The first 
three lines in Eq. 13 7|) correspond to the contribution from the quarks taking part in the condensation while the fourth 
and fifth lines corrspond to the contribution from the two light quarks with the blue color. 
Using Eq.sljSH) H38|l . the gap equations for the masses of u,d quarks are given by 

M,, - m,; = y dk^ (^°^' " '^'^ ' '''''' ^ (''°^' " ' ''''' 



(27r)= 



sin^ 6*1 - sin^ 6*!^ + |e'^ | (sin^ d_j + sin^ 0^^) - sin^ ff"^ - sin^ 6l!f^ j , (39) 



The mass gap equation for the strange quarks is given as 



Ms - m3 = / ^ ^ (1 - sin^ - sin^ Ol^). (40) 



It is worthwhile to take the zero temperature limit of Ea. l|37|l to compare with the results obtained earlier in 
Ref.0,|3|- Using the relation lim^^oo ^ ln(l + exp(— /3a;)) = —u;d{—uj), we have the contributions from the u,d quarks 
beyond chiral symmetry restored case (i.e.M„ = = Md) given as 

QudiT = 0, Af„,d = 0) = ^ - E I {VW^WT^+ VW+WT^) dk 

1=1,2 '' 

+ ^ p k'(yWT(k^W+5^)dk (41) 



where, /i± = /i± y ~ "^^^ limits of the integrations over momentum in the third integral arises from the condi- 
tions of nonnegativity of the argument of theta function for nonzero contributions while taking the zero temperature 
limit. It is reassuring that the expression is identical as in Ref.j6j|, where the two flavor colorsuperconductivity was 
investigated. 

Thus the thermodynamic potential is a function of three parameters: the two mass gaps and a superconducting 
gap which need to be minimised subjected to the conditions of electrical and color charge neutrality. The electric and 
charge neutrality constraints are given respectively as 

Qe = - - \p'' - P. = 0, (42) 

and, 

Q8 = 7fE(/^''+/'''-V') = 0. (43) 

In the above p*" = {ip^"-'' ^p'"^) = 211°- (i, a not summed) and is as given in Ea. (|23|l . Further, — X]a=i sP'"- 
It may be worthwhile to note that the neutrality conditions Ea.()42|l and Ea. H43|l can be combined to give rise to a 
simpler equation as ^Qe — -^Qs — 0, i.e. 

p''-p''+p''-pe = (44) 
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At zero temperature, for the BCS phase, the number densities of the condensing u and d quarks are the same i.e., 
p^^ = p^^, and Ea. l44() can be solved for the electric chemical potential as /i^ = — — -^/is)- K can also be shown 
by substituting this solution in Eg. 143(1 that /ig is very small compared to the chemical potentials pE and p as long 
as A//i is small 8]. 

On the other hand, for the case of gapless phase (A < one can solve Eq. (|44|l for the gap at zero temperature 

as 

--^^-{^) («) 

Ea. (|36f{lS|l and the superconducting gap equation Ea. l35() constitute the basis of the numerical calculations that we 
shall discuss in the next section. 



III. RESULTS AND DISCUSSIONS 



For numerical calculations we have taken the values of the parameters of NJL model as follows. The cutoff A 
and the scalar coupling Gg are chosen by fitting the pion decay constant f^^ — 93 MeV and the chiral condensate 
(uu)^^^ = —250 MeV—{dd)^/^. In case of non zero current quark masses, the additional parameters can be fixed from 
the values of pion and kaon masses. In the following, however, we shall choose the current quark masses for the light 
quarks as zero. This leads to the coupling constand Gs and the cut off A as = 5.0163 GeV~^ and A — 0.6533 
GeV. Similar to Ref. ,25j we take the current quark mass of strange quarks as nis = 120 MeV as a typical value giving 
"reasonable" vacuum properties. With this choice of parameters, the constituent quark masses at zero temperature 
and density are given as Mi — 0.313 GeV=M2, and for strange quark M3 = 0.541 GeV. These values are similar to 
those obtained in Ref. |2Cl(| . where the parameters have been fixed by fitting vacuum masses and decay constants of 
pesudoscalar mesons. 

Let us begin with the discussions of results without imposition of charge neutrality conditions. At zero temperature, 
the behaviour of the gap parameters as a function of baryon chemical potential (three times the quark chemical 
potential) is displayed in Fig. 1-a. We may point out here that these solutions for the gaps corresponds to the solutions 
for which the thermodynamic potential is minimised. In fact, in general, for certain values of the chemical potential 
there can be several solutions of the gap equations, but we have chosen the ones which minimise the thermodynamic 
potential. 

At low chemical potentials pB < Pi ~960 MeV, the diquark gap vanishes and the masses of the quarks stay at their 
vacuum values. The entire region below pb = Pi corresponds to vacuum solution and has zero baryon number. At 
Pb = Pi SL first order phase transition takes place and the system is a two flavor color superconductor. The diquark 
gap jumps from zero to about 80 MeV at this point. At the same point, the masses of u and d quarks drop from their 
vacuum values to zero. 

A first order transition is characterised by the existence of meta stable phases, the equivalent of e.g. oversaturated 
vapour. The masses corresponding to these metastable phases are the nontrivial (M ^ 0) solutions of the mass gap 
equation but have higher thermodynamic potential as compared to the solution corresponding to the stable phase. 

With increasing pB the superconducting gap increases until it reaches a maximum sd pB ^ 1565 MeV. The baryon 
number density becomes nonzero in the superconducting phase. At pB = Pi, jumps from zero to 0.34 fm~^ (around 
twice the nuclear matter density). The strange quark density remains zero till pB ^ 1630 MeV as strange quark 
mass remains at its vacuum value in this regime. Just beyond this point, the constituent mass of strange quark 
drops to about half its vacuum value and then decreases slowly. The density of strange quarks also becomes nonzero 
beyond pB — 1630 MeV. This will have important consequences regarding charge neutrality conditions as we shall 
discuss little later. A flavor mixing interaction probably could lead to a reduced value of the strange quark and hence 
to a lower threshold for nonzero density of strange quarks. At larger chemical potential, the results become cutoff 
dependent and the diquark condensate decreases with pB ■ We also note here that contrary to the vector- interactions 
of Ref.Q, we do not find a window in chemical potential where both chiral and the diquark condensate coexist. 

In Fig. 1-b we show the temperature dependence of masses at zero density. Chiral symmetry is restored for light 
quarks at temperature about 195 MeV. The sharp first order transition of zero temperature becomes a second order 
transition at zero baryon density as is reflected in the smooth variation of the mass which is proportional to the order 
parameter {'tljip). 

In Fig|21 we show the resulting phase diagram for the chiral transition in the temperature and quark chemical 
potential (1/3/xb) plane in the present model when charge neutrality conditions are not imposed. The middle line is 
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FIG. 1: Gap parameters when charge neutrality conditions are not imposed. Fig. 1-a shows the gaps at zero temperature as a 
function of baryon chemical potential. Fig. 1-b shows the gaps at zero baryon density as a function of temperature. Solid curve 
refers to masses of u and d quarks, the dotted curve refers to mass of strange quarks and the dashed curve corresponds to the 
superconducting gap. 



the critical line and corresponds to the state where the two nontrivial solution of the gap equations have the same 
thermodynamic potential. Along this line the thermodynamic potential has two minima of equal depth separated by 
a barrier and the barrier height decreases with temperature. At the tricritical point this barrier vanishes. Beyond this 
point there is only one solution for the gap equation and the transition is second order. This critical point {Tc, fic) 
turns out to be (74.9,285 ) MeV. The other two lines are the spinodal lines constraining regions of spinodal instability. 
We might note here that for this transition the strange quarks do not play a role as their mass is much too large to 
contribute to the dynamics. 

In FigOwe display the density dependence of the superconducting gap for various temperatures. Similar to various 
BCS-type calculations, this transition is of the second order. The critical temperature is about half of the gap value 
at zero temperature as is in ECS case. This phase transition is a second order phase transition with the transition 
temperature ~ 0.57A(r = 0). 

We next extend our discussions to the case where the charge neutrality conditions {fj,E=Q—IJ'a) are imposed. In 
Fig2]we show the dependence of masses when these conditions are imposed. For low baryon chemical potential, we 
see that the d quark masses start decreasing earlier than that of u- quark. Charge neutrality conditions force the d 
quark density to be larger (almost twice) than that of u-quark density making their masses vanishing earlier. 

At higher chemical potentials, strange quarks help in maintaining the charge neutrality condition. We might note 
here that, in general there could be several solutions of the mass gap solutions. A solution which may be free 
energetically favorable when chrage neutrality condition is not imposed can ingcneral be disfavored when charge 
neutrality condition is imposed. A typical example is shown in Fig|Sl for strange quark mass gap for electrically 
neutral normal (A — 0) u-d-s matter. 

In Fig. |S1 we display graphical solutions for strange quark mass when electrical charge neutrality condition is 
imposed. The two figures - Fig.5-a-b correspond to two different values for chemical potentials ^b- Each value 
of electron chemical potential given in x axis then gives the value of the strange quark chemical potential as Hs = 
fiB — (1/3)/xe. The dotted lines in these figures gives the solutions of the mass gap equation for the strange quark for 
each fj,s. As may be seen there can be three different solutions for given baryon chemical potential and a given value 
of electric chemical potential. We have also plotted here the charge neutrality line i.e. strange quark mass parameter 
as a function of electric chemical potential such that the total electric charge is zero. Thus, along the solid line the 
charge neutrality condition is satisfied. Along the dashed line mass gap equation is satisfied. The intersection of 
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FIG. 2: Phase diagram in the (ji, T) plane. Middle curve is the critical line and the outer lines are the lower and upper spinodals 
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FIG. 3: Superconducting gap as a function of baryon density for different temperatures. Here charge neutrality conditions 
(i.e.yE~0—fj,s) are not imposed. 



the charge neutrahty hne with the curves satisfying the gap equation defines the desired solution satisfying both gap 
equation and the charge neutrahty condition. When charge neutrahty condition is not imposed, the one which has 
minimum free energy is the required solution. However, this need not be the solution when charge neutrality condition 
is imposed. Of the two cases, shown in Fig. \E\ for /is = 1528 MeV (Fig. 5-a), the solution with masss about 520 MeV 
is the one which satisfies the charge neutrality condition. For /is = 1606 MeV, (Fig. 5-b) Mg ~ 266 MeV satisfies the 
charge neutrality condition. We might note here that, when charge neutrality condition is not imposed, the branch 
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FIG. 4: Masses of u-quarks (solid line), d-quarks (dotted line) and s-quarks (dashed line) as a function of baryon chemical 
potential when charge neutrality conditions are imposed. 
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FIG. 5: Solutions of strange quark mass gap equation and charge neutrality condition at zero temperature. In all the figures 
the electrical charge neutrality line is denoted by the solid line and solutions of the mass gap equation for the strange quark is 
denoted by dotted line for baryon chemical potential hb = 1528 MeV (Fig. 4-o.) i^n-d hb ~ 1606 MeV (Fig 4-b) 



with Ms ~ 520 MeV is free energetically preferred solution for baryon chemical potential as large as 1630 MeV as 
we mentioned earlier. The other two solutions although free energetically unstable when charge neutrality condition 
is not imposed one of them becomes the preferred solution when charge neutrality condition is imposed. Similar thing 
happens when we impose charge neutrality conditions for the superconducting gap as we shall see shortly. 

We have also verified that in presence of both electric and color charge neutrality conditions, we do not have a 
simultaneous diquark and quark antiquark condensation. 

We compute the thermodynamic potential numerically as follows. For a given temperature and quark chemical 
potential, the thermodynamic potential is minimised with respect to the strange quark constituent mass after solving 
the superconducting gap equation Ea. (|35fl selfconsistently. The values oi he and are varied so that the charge 
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FIG. 6: Superconducting gap at different temperatures as a function of the baryon chemical potential, when the color and 
electrical charge neutrality conditions are imposed. 



neutrality conditions Ea. (|42|l and Ea. H43|l are satisfied. The resulting superconducting gap is shown in FigEl for 
different temperatures. 

As may be seen for smaller values of the baryon chemical potential (/is < 1470 MeV) we have smaller values of 
the gap which reaches a maximum about 80 MeV at hb = 1390 MeV at zero temperature. Then it decreases and 
vanishes at fiB — 1470 MeV— ■ We have also verified that the difference between the value of the gap as obtained 
numerically satisfying the charge neutrality condition and the the approximate expression given in Eg. 145(1 equation 
is less than a percent for the gapless phase. The gap remains zero till /is = /^a = 1530 MeV. In fact, in the region 
between fi2 and ^3 we do not have any real solution for A. At /ia the system jumps from the normal phase to the 
BCS phase with a gap about 132 MeV. The baryon number density jumps from pB — O.OGfm""^ to 1.54 fm~'^ i.e. from 
6 to 9.5 times the nuclear matter density. The interval in the baryon chemical potential between the gapless and BCS 
phase within which it is normal quark matter decreases with increase in temperature and then disappears as may 
be noted in Fig.fEl The sharp transition between gapless phase to the BCS phase as a function of baryon chemical 
potential at smaller temperature becomes a smooth transition as temperature increases as may be seen from Fig. 6-b. 
We also would like to note that the starnge quark number density at zero temperature is zero in the gapless phase. 
This is in contrary to observation made in Ref.0 where a vector vector point interaction was considered. However, 
as temperature is increased, the number densitiy of the strange quarks become nonzero. E.g., at a temperature of 50 
MeV, strange quarks contribute upto about 15 % of that of density of u-quarks. 

Let us next discuss now some of the characteristics of the gapless modes which occur for smaller values of the baryon 
chemical potential {fj,B < M2)- In fact, in this region, in general, there are two solutions to gap equations for given p. 
and S^. One is the usual BCS solution with a larger value of the gap and the other one with a smaller value of the 
gap. However, it so happens that in this range of the chemical potential, the solution of the gap equation with larger 
value of the gap (the BCS solution) does not satisfy the charge neutrality condition. The situation here is similar to 
the case for strange quark mass when charge neutrality condition is imposed. 

At zero temperature, gapless modes occur when the gap is less than half the difference of the chemical potential 
6fi of the two condensing quarks. It is easy to show also that the excitation energy uj2 of the d-quark vanishes at 
momenta — p. — ^ — and = /2 + \J — A^ at zero temperature. In this phase the number densities of 
the condensing u and d quarks are not the same. Let us note that the number densities of u-red quarks are given as 



= (^" V") = / rfk(f " - Fr") = 9 



(46) 
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FIG. 7: Occupation number as a function of momentum. Solid line corresponds to u quarks and dotted line corresponds to d 
quarks 



where 



for the quarks and 



F" = sm^ 9l-Ul-^] {1- sin2 6l - sin^ 9l) 



" - sin^ e\ 



2 Ql 



,2 q2 • 



(47) 



(48) 



for the antiquarks. Interchanging indices 1 and 2 we shall have the expression for the number density of the d 
quarks taking part in condensation. At zero temperature, the distribution function sin^ 9\_ vanishes for the u quarks 
as coresponding energy uoi — ui — 5^j, in always positive (let us remember that is negative at zero temperature) 
where as the distribution function for d quarks is nonzero between /i_ and In fact, the number densities at zero 
temperature are given as 



(27r)3 



2 



~^2)) - ( 1 - ^ 



(49) 



for the u quarks. Similarly the density of the d quarks taking part in condensation is given at zero temperature by 
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FIG. 8: Number densities of u quarks (solid) and d quarks (dot-dashed) participating in superconducting phase at zero (Fig 8. a) 
and at finite temperatures (fig. 8.b). 
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FIG. 9: Excitation energies as a function of chemical potential. 

The integrands of Eq.s (|49() and (|50|l which are the occupation numbers in the momentum space are plotted in Fig|7| 
It is easy to see from the expressions in Ea.s (|49|l and l|5U|l that except in the interval in the momentum 

space, the theta function does not contribute and the distribution is like the BCS distribution. The occupation 
number distribution is a 'dented' theta function with the 'dent' ~ A around the average fermi surface. In the region 
between /i_ and /^+, the occupation numbers resemble like that of normal matter as pu becomes almost zero but for 
the vanishingly small antiparticle contribution. On the other hand, the integrand of pd becomes unity in the same 
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limit for the antiparticlc contributions. This is shown in FigEI It foUows therefore that the total number of u-quarks 
and d-quarks (integrated over the momentum) are not the same. The difference in their number densities is in fact 
given by 

5psc = P2 - Pi = J dk [(sin^ 0l - sin^ 0l) + (sin^ 0l - sin^ 0^)] . (51) 

Thus at finite temperatures the Spsc will always be nonzero as long as the difference in the chemical potential is 
nonzero. At zero temperature this difference will be nonzero only for the gapless solution i.e. for gaps less than half 
the magnitude of the difference in chemical potentials. This is infact given by 5psc{T = 0) = 2/(37r^)(/ii^ — p!t.). In 
Fig|Sl we plot the number densities of u and d quarks taking part in the condensation. At zero temperature, the 
number densities become the same in the BCS phase. However, because of thermal distributions, the two densities 
do not become exactly identical even in the BCS phase at finite temperatures. 

In Figl^we plot the dispersion relations - the excitation energies of the quasi particles as a function of momentum 
for quark chemical potential pq — (l/3)/is ~ 379MeV and at zero temperature. Let us note that the dispersion 
relations are given by wi = uj — 5^ and uj2 — uj + 5^. The gap in this case turns out to be A = 59 MeV and the 
difference of the chemical potential turns out to be (5^t ~ —66 MeV at zero temperature. Far from the pairing region 
Jl ~ 354 MeV, the spectrum looks like that of usual BCS type dispersion relations. Of the two excitation energies, 
uji shows a minimum at the average fermi momentum p with a value a;i(|k| = fi) — A + (5^ ~ 125MeV. On the 
other hand, UJ2 becomes gapless {uj2 — 0)at momenta /i_ and In this 'breached' pairing region |26| we have only 
unpaired down quarks and no up quarks. 

In Fig^lwe show the temperature dependence of the gapless modes. We have taken the quark chemical potential 
as 379 MeV as earlier. Initially the gap increases and then decreases monotonically and vanishes at Tc — 48MeV i.e. 
Tc = 0.81A(T = 0). Clearly, this shows a departure from the BCS behaviour with Tc/A{T = 0) ^ 0.567. Further, it 
is also seen that for temperatures between 28 MeV to 32 MeV, the gap is larger than difference of chemical potential 
and UJ2 no longer is exactly gapless in between these temperatures. The excitation energies however are much smaller 
than temperature and are almost similar to the excitation energies between /i_ and at zero temperature. For 
temperatures beyond 32 MeV, the gap becomes smaller than and we have again the gapless modes occuring at 
corresponding /i_ and 
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FIG. 10: Temperature dependence of gap parameter in the gapless phase (solid line). \S^\ as a function of temperature is shown 
as the dotted line. 

It may be worthwhile to mention here that gapless modes in superconductivity were known theoretically in the 
context of superconducting matter with finite momentum ,23J as well as in condensed matter systems with magnetic 
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impurities ;27i . Recently it has been investigated for cold fcrmionic atoms !26L |23 | . Gapless modes in the context of 
quark matter has been first proposed in Ref.|2J| for color flavor locked matter. However, this corresponded to a meta 
stable phase. Gapless modes in neutral quark matter were first emphasized for the two flavor colorsuperconductivity 
in Ref. and for the 2SC+S quark matter in Ref. Stable^apless modes for color flavor locked matter has 

been first proposed in Ref. '29^ and has been confirmed in Ref. 30] in a more general structure for the gap. The 
temperature dependence of the CFL gapless modes has been studied in Ref. j3lj . 

As the chemical potential is increased beyond 1530 MeV at zero temperature the solution for the gap jumps 
to a higher value of about 130 MeV almost similar to the case when charge neutrality condition is not imposed. 
This corresponds to the usual BCS solution. In this case, the numbers of u quarks and d quarks participating in 
condensation are the same. The charge neutrality condition however is maintained by the third color, the electron as 
well as the strange quarks. One essential effect of including strange quarks is that the electron density starts decreasing 
for chemical potentials greater than the strange quark mass as the strange quark can carry the negative charge to 
maintain electric charge neutrality condition. This has the effect that the lowest excitation energy UJ2 ^ + in 
the BCS pairing case becomes large due to both the large value of the gap as well as the smaller magnitude of the 
electron chemical potential. 

IV. SUMMARY 

We have analysed here in N JL model the structure of vacuum in terms of quark-antiquark as well as diquark pairing 
at finite temperature. The methodology uses an explicit variational construct of the trial state. The gap function 
as well as the thermal distribution functions are also determined variationally. To consider neutron star matter, we 
have imposed the color and electric charge neutrality conditions through the introduction of appropriate chemical 
potentials. It has been noted and emphasized earlier that projecting out the color singlet state from color neutral 
state costs negligible free energy for large enough chunk of color neutral matter Q . 

We observe that the mass of the strange quark plays a sensitive role in maintaining charge neutrality condition. 
Further, we observed that solutions for the gap quations which may not be free energtically preferble solutions, can 
be the preferred solutions when constraints of charge neutrality conditions are imposed. We also find that with 
increase in chemical potential, the quark matter has a transition from a gapless phase to the BCS phase through 
an intermediate normal quark matter phase. This interval between the gapless and BCS phase decreases and then 
disappears with increase in the temperature. We also find that there are no strange quarks in the gapless phase at 
zero temperaure contrary to the observation in the vectorial point interaction model as in Ref. 7]. However, at higher 
temperatures only the contribution of strange quarks becomes significant even in the gapless mode. We have looked 
into the gapless to BCS transition as a function of density for various temperatures. The sharp transition at zero 
temperature becomes a smooth transition at higher temperatures as density is increased. 

We have focussed our attention here to the two flavor superconducting phase with strange quarks. The variational 
method adopted can be directly generalised to include color flavor locked phase and one can then make a free energy 
comparison regarding the possibility of which phase would be thermodynamicaly favourable at what density. This 
will be particularly interesting for cooling of neutron stars with a CFL core. We have considered here homogeneous 
phase of matter. However, we can also consider mixed phases of matter with matter being neutral on the average. 
Surface tension between the sucrconducting phase and the normal quark phase will become an important factor in 
determining the stability of this mixed phase |32i ■ Some of these problems are being investigated and will be reported 
elsewhere l33l. 
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